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Abstract 

We prove Knudsen's law for a gas of particles bouncing freely in a two dimensional pipeline 
with serrated walls consisting of irrational triangles. Dynamics are randomly perturbed and 
the corresponding random map studied under a skew-type deterministic representation which 
is shown to be ergodic and exact. 
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1 Introduction 

In his nowadays classical studies on the kinetic theory of gases, the Danish physicist M. Knudsen 
experimentally observed that, no matter how an inert gas was injected into a pipeline, the direction 
in which a molecule rebounds from the pipeline's solid wall is asymptotically independent of its 
initial trajectory. That is, the fraction of particles leaving the surface in a given direction is propor- 
tional to cos (8), where 8 is the angle that such a particle's trajectory defines, measured with respect 
to the normal to the surface ([H]). This behaviour is referred to as the Knudsen's (cosine) law 
ever since. In the experiment, the gas is injected at a very low pressure so that interactions between 
particles are negligible. 

The physical justification of Knudsen's law is the following. First of all, one assumes that 
particles bounce at the pipeline's wall elastically. This means that the energy of a particle is 
preserved in a collision which, in turn, implies the refection law: the angle that the incident direction 
forms with the normal to a surface, 8i n , equals the angle 8 re f of the reflected direction, 8i n = 8 re f. It 
has been proved that this law is valid in a first approximation if we do not take into account thermal 
effects. Since we also assume that particles do not interact with one another, we must conclude 
that the microscopic irregularities on the pipeline's surface are responsible for the destruction of 
any particular pattern in the original gas distribution. Indeed, even if we assume that the bounces 
are perfectly elastic, microscopic holes in the boundary of the pipeline and imperfections in relief 
are dimensionally comparable to the molecules of the gas and have therefore a disruptive (i.e., 
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unpredictable) effect on particle collisions: after many bounces, we are in the so-called Knudsen's 
regime in which the reflected direction is independent of the incident one. However, this argument 
does not explain by itself why the reflected angles are distributed according to the cosine law. The 
theory of billiards helps clarify this point ([H], [7)). 

The irregularities of the pipeline surface can be reasonably modelled as cavities or microscopic 
cells with a dispersive geometry such that, once a particle has entered one of them, it comes out 
of it with a rather arbitrary direction, even if all the collisions inside that cavity are elastic. If 
the pipeline is made of a uniform material, we can model one such cell as a billiard table and the 
pipeline wall as an infinite row of such billiard tables that a particle moving freely enters and exits. 
With this description, Knudsen's cosine law can be then seen as a consequence of the fact that, for 
sufficiently dispersive billiard geometries, the Liouville measure is the unique measure preserved by 
the billiard flow (see [3] for more details). Unfortunately, this need not be the case for polygonal 



(i.e., non-dispersive) tables such as Fig. 2.1 So polygonal geometries require a slightly different 
approach. 

Since the characteristic size of this cell is infinitely small compared with the diameter of the 



pipeline, the exact position along the open side of the cell (dotted line in Fig. 2.1| at which 
the particle enters a cell coming from the previous one is considered to be randomly distributed 
with uniform probability. Following jB], we will refer to these billiards as random billiards. More 
concretely, the dynamics of a particle bouncing inside the pipeline are given by the first return 
map to the open side of a the billiard cell, which defines a Markov process. Thus, the dynamics 
are characterised by a transition operator K (9, A) such that, for any reflected direction 9 £ [0,7r], 
K(9,A) is the probability that a molecule takes a direction in A C [0,7r] after the next rebound. 
These concepts will be reviewed in Sections [3] and [4] In this formalism, Knudsen's law can be 
mathematically written as follows. 

Let v be the initial angle distribution with which the gas is injected into the pipeline. That 
is, if A C [0, 7r], v(A) represents the proportion of particles that will hit the pipeline's wall for the 
first time with an incident angle 9 £ A. Denote by z/™) the distribution after n collisions. In this 
context, Knudsen's cosine law can be expressed in two different ways. The strong Knudsen's law 
claims that, for almost any initial distribution v, after many collisions v^ n > converges to the cosine 
distribution, 

V W(A) — > 1 / cos{9)d9, AC[0,tt]. (1.1) 



A 



On the other hand, the weak Knudsen's law states that if, after any collision, we count the 
number of particles reflected with a given angle 9 and then divide by the total number of particles, 
this quantity is proportional to cos((9). Explicitly, 



n 

-5>W(A) -> - / cos(9)d9. (1.2) 

i—1 



A 



Obviously (1.1 1 implies (1.2 1, so the strong Knudsen's law implies the weak one. 

Deterministic dynamical systems can be sometimes understood as idealizations of real systems, 
the latter usually subjected to negligible perturbations in a first approximation. The use of a 
dynamical approach to study a gas is one such example. In this paper, we will study the properties 
of the random billiard obtained as a result of modelling the rough surface of the pipeline by a 
zig-zag geometry as shown in Fig. |2.1| Our billiard cell is an isosceles triangle with open side pq 
defined by an angle a. This model was introduced in [7J but has not been studied in depth yet. For 
example, in [7] nothing is said about the ergodic properties of this billiard table when a/2ir ^ Q (it 
is shown that it is not ergodic under the Liouville measure if a/2n £ Q though). Thus it is not at 
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all clear if the dynamics resulting from a random billiard as in Fig. 2.2 will in fact follow Knudsen's 
law (weak or strong). In this paper, we will prove that it is actually the case. More concretely, the 
main contributions of our paper are the following: 

1. Little is known about the ergodicity of the flow of deterministic billiards given by irrational 
polygon tables Remarkably, we show that, for a concrete example, randomising the 
billiard in a sensible way (i.e., choosing the point at which the particle enters the billiard ran- 
domly but uniformly distributed on one side of the table) implies its ergodicity and exactness 
with respect to the Liouvillc measure. 

2. Unlike other attempts to prove (weak) Knudsen's law, where one usually assumes very irreg- 
ular geometries at a microscopic level responsible for the dispersive effects, we deal with an 
extremely simple pattern. This has additional advantages as a simpler model capturing the 
main features of a system can be simulated and developed more easily. 

3. In the literature, ergodicity of random billiards relies on the ergodicity of the first return map 
(see [7]). We use here a different approach and take advantage of a skew- type representation 
of random maps recently introduced in [T] . The techniques used to prove the exactness of our 
model can be certainly adapted to study other (random) billiard tables and random maps in 
general. 

4. Finally, we prove that the strong Knudsen's law holds for our model. As we mentioned before, 
even if one assumes dispersive billiard geometries, one will only prove that the first return 
map is ergodic, i.e., the weak law. 

The paper is structured as follows: in Section [2] we introduce the random billiard behind our 
model. We recall some basic definitions and properties of billiards in Section [3] and the concept 
of random billiard in Section |4j In Section [5j we review the skew-type representation of random 
maps introduced in [T] and show the relationship between this representation and the dynamical 
evolution of absolutely continuous measures. We prove that this skew-type representation is exact 
in Section [7J which uses some auxiliary results presented separately in Section [6] Finally, in Section 
[8j we illustrate our results with numerical simulations. 

2 The Model 

Suppose we have a gas of non-interacting particles moving freely inside a pipeline. Our model 
consists of a two dimensional infinite pipeline whose rough walls are modelled as a sequence of cells 



built as a juxtaposition of a fundamental cell (Fig. 2.2). In this section, we are going to fix the 
geometry of such cell and introduce the main hypothesis behind the dynamics of the particles. The 
content of this section is extracted from [TJ. 

Both the geometry and dynamics of our model are summarised in Figure [2~T] We assume that 
the walls of the pipeline are not smooth but describe a serrated regular pattern built from a fixed 
isosceles triangle with an open side. We characterise that triangle by one of its angles a. The 
dynamics of a particle moving freely inside the pipeline are as follows: the particle enters one of 
the cells with some angle 9 in and bounces on its sides elastically, that is, according to the reflection 
law. The particle eventually comes out of the cell with a given direction 9 outl then crosses the 
pipeline, and reaches another cell on the opposite wall, and the process recurs. Since the diameter 
of the pipeline is several orders of magnitude bigger than the characteristic length pq of the cell, it 
is plausible to think that, every time the particle enters a new cell through the open side, it does 
so at a point uniformly distributed on pq. 
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Figure 2.1: Particle bouncing in a pipeline with serrated triangular boundaries. 

In the literature, one encounters examples of particle dynamics on bounded domains where the 
random perturbation is introduced in the observations of 9 out every time the particles bounces off 
the wall (see for example [5]), which implies that collisions are no longer elastic. This perturbation 
is explained as a consequence of the roughness of the wall. Observe that, in our model, randomness 
is introduced in a completely different way. We assume clastic collisions but, instead of considering 
dispersive geometries responsible for the random bounces, the noise is introduced in a sensible way 
without modifying the (rather elementary) geometry of the problem. 

The trajectory that a particle describes after it comes out of a billiard cell until it enters another 
one on the opposite wall is completely irrelevant for dynamical purposes. In practice, we can 
better understand our dynamical system as the closed billiard table obtained by fixing together two 
triangular cells so that they share the same open side as in Figure [iO| When the particle crosses 
that open side, it retains its direction but it is assigned a different position on pq according to a 
uniform law. 




Figure 2.2: Fundamental cell. A particle crosses the open side pq at a uniformly distributed random 
point. 

The incoming angle 9i n (resp. outgoing angle) is then the angle in [— tt/2, tt/2] with which 
the particle enters (resp. leaves) the cell measured with respect to the normal to the open side pq. 
One often shifts [— 7r/2,7r/2] by tt/2 and measure 8 in [0, ir], i.e. with respect to horizontal. In [7], 
for any incoming angle e [0, 7r], Feres gives the different possible directions with which a particle 
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may come out of the fundamental cell (Fig. 2.2 1. They are four, given by maps t,; : [0, tt] — » K 
defined as 

Ti (0) = 6> + 2a, t 2 (0) = -0 + 2tt - 4a 

t 3 (0) = - 2a, 74 (0) = -0 + 4a, (Ml) 



and any 7j applies with certain probability p t . The probabilities associated to the maps (Ml) are 
as follows: 



where 



Pi (0) 



P2 (0) = 



P3 (0) = 



P4(0) 



1 0G[O,a) 

u a (0) e [a, rr - 3a) 

2 cos(2a)u 2a (0) e [rr - 3a, tt - 2a) 
G [tt - 2a, tt] 

0e[O,7r-3a) 

u a (0) - 2 cos(2a)u 2a (0) G [tt - 3a, tt - 2a) 

w a (0) G [tt — 2a, 7r — a) 

G [7T - a, tt] 

G [0, 2a) 
2cos(2a)u 2a (-0) 0e[2a,3a) 
u a ( — 0) G [3a, 7r — a) 

1 G [7r — a, 7r] 

G [0, a) 

u a {-ff) 0e[a,2a) 

u a (-0) - 2cos(2a)u 2a (-0) G [2a, 3a) 

0G[3a,7r], 



u a (0) = s 1 + 



tana 
tan0 



(M2) 



(2.1) 



and a is assumed to be smaller than 7r/6. As we will see later, this family of maps and probabilities 
is all we need to prove Knudsen's law. 



3 Preliminaries on Billiard Dynamics 

This section aims at recalling the main concepts of the theory of (deterministic) billiards and one 
of its key results, namely, that the billiard map preserves the Liouville measure. The content of 
this section was extracted from [3], which the reader is referred to for an exhaustive exposition on 
billiards. 

A billiard table D is the closure of a bounded open domain D C K 2 whose boundary dD = 
Ti U • • • U r m is a finite union of smooth compact curves (C l , I > 3). Ti, T n are called the walls 
of the billiard table and are assumed to intersect each other only at their endpoints or corners. 
For any i = 1, ...,n, is defined by a C l map fi : [aj,6j] C K —> M 2 such that the second 
derivative f" either vanishes or is identically zero. This condition prevents the wall from having 
inflection points or line segments. A wall such that /" ^ is called focusing or dispersing if /" 
points inwards or outwards D, respectively. Otherwise /" = and the wall is called flat. In this 
framework, we consider the dynamics of a point particle moving freely within D. 

Let q G D denote the position of a free particle moving within a billiard table. We suppose that 
the particle bounces at dD elastically. That is, if it collides with the wall T at a point p that is 
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not a corner, the incident angle that the velocity v forms with the normal vector n to T at p 
equals the angle of reflection 9 re f between n and v after the collision, i.e. 9- m = 9 re f. This implies 
that v has constant norm, which we will assume equal to 1 for the sake of simplicity. Under these 
assumptions, the billiard flow $ is the flow that this particle defines on the phase space D x S 1 . 
Since the particle follows a straight trajectory between two consecutive collisions, the dynamical 
information of the system is contained in the geometry of the boundary dD and how the particle 
bounces off it. Therefore, instead of studying the billiard flow, first one introduces the cross-section 
M = dD x [— 7r/2,7r/2] as the set of all postcollisional velocity vectors, where 9 € [— 7r/2,7r/2] 
measures the angle between v and the normal vector n pointing to the interior of D, and then 
considers the first return map F that the billiard flow induces on M. That is, given a point of M, 
F gives the next position at where the particle collides and its velocity after that collision. The 
map F is often called the billiard map. 

It is a well known result in the theory of billiards that the billiard map preserves the Liouville 
measure A<8>/z on M, where A is the Lebesgue measure on dD and /i is the measure on [— 7r/2, 7r/2] 
given by u(A) — | J A cos(9)d9, A C B ([— tt/2, 7r/2]) ([21 Lemma 2.35]). However, except for a few 
general results, little is known about the ergodicity of billiard maps with respect to the Liouville 
measure for explicit examples of boundary geometries. For example, dispersive billiards are ergodic 
(0 Theorem 6.20]) while regular polygons are not. 

Finally, as we pointed out in Section[2j it is customary in the literature to translate [— tt/2, tt/2] 
and measure 9 in [0,7r], in which case \x is given by fi(A) — \ J A sin(9)d9, A C B ([0, tt]). We will 
follow this convention throughout the paper. 

4 First Return Map. Random Billiards 

In this section we are going to get back to the concept of random billiard and its different repre- 
sentations. 

Let D be a billiard table with boundary dD = T\ U • • • U T m . Remove one of the walls of dD, 
for example T±, so that we obtain a billiard table with an open side. Without loss of generality, 
we can assume that the open side is just a line segment pq as in Figure |2.2| This identification 
can be carried out by the map /i : [01,61] C K — > M. 2 defining Ti. On the other hand, suppose a 
particle enters the billiard from the open side at x e pq with an angle 9i n € [0, tt] with respect to 
Ti. Denote by 1 $ x (9in) G [0,tt] the angle associated to the velocity of that particle when it returns 
to the open side for the first time. The first return map (to the open side) is then the map 
T : pq x [0, tt] —> pq x [0, tt] induced from the billiard map F such that T (x, 9 in ) = (y, 9 x (9i n )), i.e., 
the particle leaves the billiard at y 6 pq with an angle 9 out = ^ x {^in)- 

The Liouville measure A ® /i on dD x [0, tt] induces a measure onpgx [0, tt] by restriction in a 
natural way. We will continue referring to this measure as the Liouville measure and denoting it by 
A ® (J,, where A now stands for the Lebesgue measure on pq. It can be proved that the first return 
map T also preserves the Liouville measure ([U Exercise 1.5.1]). Furthermore, 

Lemma 1 If the billiard map F : dD x [0, tt] — > dD x [0, tt] is ergodic with respect to A ® fJ>, so is 
the first return map T : pq x [0, tt] — > pq x [0, tt] . 

Proof. Let A € B(pq) ® B([0, tt]) be an invariant set, i.e., T -1 (A) = A. This set can be naturally 
regarded as a measurable set of dD x [0, tt] since pq C dD is identified with the wall Ti. Since F is 
ergodic, U n >i-^' _ ™ (^) nas ^ uu ^ ® ^la_D X [o ^-measure, which means that 

A ® V\ 9DX p.,] ((Un>iF- n (A)) P| (pq x [0, tt])) = 1. 
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If y € F n (A) f] (pq x [0,7r]), then y = F n (x) for some x 6 A and, from the very definition of 
the first return map that, y = T~ mn (x) for some m„ < n. Therefore, 

u F ~ n ( a )) n >< to, *}) = u >1 T ~ m ( a ) = a > 

^-^ n>l / 1 1 ^-^ m>l 

where the last equality follows from the invariance of A. Consequently, A <8> Mlpq X [o w ] (^) = 1 an d 
T is ergodic. ■ 

Let (pq, B(pq), A) be the probability space built upon pq with the normalised Lebesgue measure 
A on the Borel er-algebra B(pq). The random billiard associated to D is then the time-discrete 
random dynamical system {9 X : [0, tt] — > [0, tt] \ x £ pq}, where pq is regarded as a probability space. 
Roughly speaking, the random billiard thus built models the outgoing state a free particle that 
entered the billiard through the open side at a point uniformly distributed along pq. When regarded 
as a probability space, pq will be denoted by SI. 

The random billiard {^ x : [0,7r] — > [0,7r] | x S SI} defines a transition probability kernel K : 
[0, tt] x B ([0, tt]) -> [0, 1] given by 

= A ({x e SI | ^(0) e A}), where z € [0,tt] and A e B([0,tt]) . (4.1) 

We say that v is invariant with respect to : [0, 7r] — > [0, 7r] | x G SI} if 



i/(A) = / K{6, A)dv{0) for any Ae B ([0, tt]) 

and that the random billiard is ergodic with respect to v if T is ergodic with respect to A <g) v. 
For example, /i(A) = 5/4 sin#cZ0, ^4 G i3 ([0, ir\), is an invariant measure because A <E) /i is invariant 
with respect to the first return map T ( 7, Proposition 2.1]). Moreover, since dispersive billiards 
are ergodic with respect to the Liouville measure, the corresponding random billiard is also ergodic 
with respect to fi (Lemma [l}. In this situation, Birkhoff's Ergodic Theorem implies that 

-Vi a (tt 2 oT") — ► I f sm(6)d9 a.s., (4.2) 

i=l JA 

where 7T2 : pq x [0, tt] — > [0, tt] is the projection onto the second factor. That is, the average number 
of particles reflected with some angle 8 £ A converges to \ J A s'm(9)d9 so (weak) Knudsen's law 
holds for dispersive billiards. 

In general, for non-dispersive billiards, little is know about the ergodicity of the first return 
map. This is indeed the case of the random billiard introduced in Section [2] (Fig. 2.2 1. Recall 



that, as far as we know, determining whether a general irrational triangular billiard is ergodic is 
still an open problem ([9l Section 7]) and only a few particular examples were proved to be ergodic 
([14 ). Therefore, at this point, it is not clear at all if Knudsen's law holds for the pipeline model 
introduced in Section [2] To prove it, the first return map must be replaced with a more convenient 
representation. 



5 Skew-Type Representation of Random Maps 

Given a random billiard, the probability kernel ( |4.1[ ) contains all the dynamical information of the 
system and determines its properties. This implies that, as far as the dynamics is concerned, the 
underlying probability space SI (which in Section [4] was equal to (pq, B(pq), A)) plays a secondary 



Lazaro-Cami, Dingle, and Lamb: Knudsen's law and random billiards in irrational triangles 



8 



role. In particular, other probability spaces (fi, J 7 , P) giving rise to the same probability kernel 
are available. Indeed, if K : X x B — > [0, 1] is a transition probability kernel on a measurable 
space (X, B) , Kolmogorov's Existence Theorem on Markov processes guarantees that there exists 
a probability P defined on the Borel cr-algebra of the topological space — YiTLi-Xi, '■— X for 
all i > 1, such that the chain ^ := 7Tj is Markovian and has transition probability kernel K ([31 
Theorem 2.11]), where TTi : \^° =1 Xi — > X is the projection onto the ith factor. In this statement, X 
has to be a cr-compact Hausdorff space and B the Borel cr-algebra. Unfortunately, Il^i^i ^ s n °t 
very manageable to work with explicitly. For instance, the representation T : pqx [0, ir] — > pq x [0, 7r] 
given by the first return map introduced in Section [4] only involves finite dimensional spaces and 
hence seems more convenient. We are now going to consider another skew- type representation for 
random maps recently introduced in pQ . As we will show in Section [7j this representation will be 
crucial to prove the asymptotic properties of the random map defined by our random billiard Fig. 

221 

Let (X, B, v) be a general measure space and let ([0, 1), B ([0, 1)) , A) be the unit interval regarded 
as a probability space, where A stands for the Lebesgue measure. Later on, we will apply the results 
of this section to X = [0, tt]. For any k = 1, N, let Tfe : X X and pk ■ X — >• [0, 1] be measurable 
mappings such that {pk}k=i,...,N is a measurable partition of the unity, i.e., J2k=iPk( x ) = 1 f° r 
any x € X (compare to (M2) and (Mil in Section|2|. We define the random dynamical system 



t '. X — )■ X such that r (x) = Ti (x) with probability pi (x). The transition probability kernel of r 
is given by 



N 

K(x,A)=^ Pi (x)l A (Ti(x)). (5.1) 



This probability kernel defines the evolution of an initial distribution (probability measure) v on 
(X, B) under the random map r iteratively as 

v^:=v, v {n+1 \A) = f K(6,A)dvW(6), n > 1, (5.2) 

where A e B. A measure v is called invariant if ia 1 ) = v. 

Following [T), consider now 51 = [0, 1) x X and set Jk ■= {{y,x) € Q \ J2i<kPi( x ) — V < ~ 
^2i<kPi( x )} sucn that { Jk}k=i....,N is a finite measurable partition of f2. We define the skew-type 
representation of the random map r as the map S:[0,l)xl4[0,l)xl such that 

S(y,x) = (<p k (y, x),T k {x)) for (y,x)eJ k (5.3) 

where 

The map iff. is well defined because if (y, x) € Jk then pk(x) > necessarily. Moreover, S thus 
defined is B([0,1)) (8) Immeasurable ([TJ Lemma 3.1]) and it is a skew-type representation of the 
random dynamical system r. 

As we will see, the map S is extremely useful to study the properties of r. For example, 
v is an invariant measure of t on (X, B) if and only if A €D v is an invariant measure of S 
on ([0, 1) x X, B ([0, 1)) <g> B, A ® v) (pQ Lemma 3.2]). Moreover, sinc e w e know that the sine 



n(A) = | J A sm(9)d8 is invariant by the transition probability kernel (5.1) of our billiard table 



we conclude that fj, is invariant by r built from (M2 ) and (Ml ), and therefore A (g> [i is invariant by 



the corresponding map S. In an abuse of terminology, we will continue calling A <£> /i the Liouville 
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measure when referring to S. In this context, we will say that the random map t with initial 
distribution v is ergodic (rcsp. mixing, exact) if S is ergodic (resp. mixing, exact) with respect 
to A ® v. In Section [7j we are going to prove that the S associated to (M2) and (Ml) is exact. 



V 4 / 






h 


Jl 






w 



Figure 5.1: The sets Jj, i — 1, ...,4, for the random dynamical system defined in Section [2j 



Unlike the first return map T (Section |4j), the map S has no dynamical interpretation and is 
a purely auxiliary tool to represent and study r. Therefore, we need to establish the relationship 



between S and the transition probability kernel (5.1), which carries all the dynamical information 



This is the content of Theorem [4| Before stating that relationship, we need to introduce some 
notation and an auxiliary lemma, whose proof is included separately in the Appendix for the sake 
of a clearer exposition. 



Definition 2 We define iteratively the sets Ji 1 ...i n , where ij £ {1, N} for any j = 1, 

J n ... ln := {oj£j tn I S(u)£j in _ 1 ,S 2 (u)£j in _ 2 ,....,S n - 1 (u)£j il }. 
Lemma 3 Let I x := [0, 1) x {x} the fibre through x £ X . If I x f] Ji x ..-i n 7^ 0, then 

Theorem 4 Let v be a measure on {X, B) and A £ B. Then, 

E A ®, [1a (tt 2 o S n )} := / 1 A (tt 2 o S n ) d(\®v) = v {n) (A) . 



..,n, as 
(5.4) 



'[0,l)xX 

Proof. First of all, it is not difficult to check that 



„<•» (A) . £ /*(.)„. WPi ._, („„ (.)) ... Pn o ... o U („, o ... o T( „ W , . 

ii,... ,i n 

This can be proved iteratively rewriting i^T in terms of (auxiliary) Dirac deltas, K(x,dy) = 
E^ift («) « (V - n (x)) dy where / / (y) S (y - a) dy = / (a), f £ C (E). 
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On the other hand, since [0,1) x X = [J i . Ji 1 ...i n is the disjoint union of the sets Ji x ...i n , 

ij e {i,...,N}, 

[ l A (l 2 oS")d(A®(/)= V f l A (TT 2 oS n )d(\®v) 

J[0,l)xX i u ...,i n J Jii -in 



^2 / Ia^O ...OT in )d(\®L>) 

/ lj n (T I1 o...or,Jd(A®!/), 

„• J[0,l)xX 



Hi — j»n 



where, in the second line, we have used that 1a ("2 o S*™) = 1a (r^ o ... o n n ) on the set Ji 1 ...i n - If 
now we apply Fubini's theorem, 



Iji! i„ 1 A(T il o... or» n )d(A®i/) = / (Tii ° — ° n n (x)) dv (x) / lj. in (y,x)d\(y) 

[0,1) X.X ix •'[O.l) 

but by Lemma [3j for any fixed a; G X, 

/ in {y,x)d\{y) =p in (x) ■ ■ -Pi^T^ o • • • o n n (x)). 

J [0,1) 

Therefore, 

(t; 2 o ... o r in (x)) 1a (t 4i o ... o Tj n (a:)) cfe (x) . 



/ lA(7r 2 o5 n )d(A(8)i/) = V / Ptn {x)--- Pll (i 

i[0,l)xX i 1 ,...,i n -> x 



Corollary 5 Lei J7 S : L p ([0, 1) X AC, A ® //) -> L p ([0, l)xX,A«(j), p > 1, be the Koopman 
operator associated to S, i.e., U$f = f ° S. Let 7r 2 : [0, 1) x X — > X be the projection onto the 
second factor. If v is a measure on (X,B), then 

i/**>(A) = / U% (K-i( A) ) d(A®i/), 4e8, n e N. 

J[0,1)x[0,tt] V 2 V 1 

Let S* be the skew-type representation of the random billiard introduced in Section [2] and let 
v <C /i be a probability measure absolutely continuous with respect to the Liouville measure [i with 
Radon-Nikodym derivative f <E L 1 ([0,7r],/i). One can easily check that = n 2(f), where 

7r|(/) (w) := / (7T2 (w)) and w G [0, 1) x [0, 7r]. If 5 is mixing (or exact) then, from Corollary [5] 

^(/^(VtA))^^)^ J n* 2 (f)d(\®v) J l^ {A) d(\®fi)=p(A) (5.5) 

( |12[ Proposition 4.4.1 (b)]), which implies the strong Knudsen's law. We will prove that S is exact 
in Section [7] and give more details about the strong law in Section [8] The proof uses the fact 
that the pull-back nZClc) of a characteristic function 1q, C <E £>([0,7r]), cannot be invariant by 
S 2 := S o S (Section^). It is worth observing that the strong law, unlike the standard approach to 
random billiards available in the literature from the first return map (see |4.2[ ) , is a consequence of 
the asymptotic properties of the skew-type representation S. 
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6 Properties of S z 



Let S be the skew- type representation of the random map (Ml) and (M2). In this section, we are 
going to show that S 2 = SoS : [0, 1) x [0, tt] — > [0, 1) x [0, tt] cannot leave invariant any characteristic 
function tt^ (lc), where C £ B ([0, n]) has probability < \i (C) < 1 and tt 2 : [0,l)xX->I denotes 
the projection onto the second factor. This result will be used in the next section to prove that S 
is exact. First, we will show that our model exhibits a very useful symmetry. 



Proposition 6 Let S be the skew-type representation of the random map (Ml) and (M2). Let 
g : [0, tt] — >• R be a function such that ^(g) is S 2 -invariant, i.e., 



^2(9) = ^2(9) s 2 , 

and let <f> : [0, 7r] — > [0, tt] be defined by 9 n- tt — 9 . Then ^(ff 



(6.1) 



is also S -invariant. 



Proof. It is not difficult to realise from ( M2 1 that 



pi o cj) = p 3 , p 3 o cf) = j3 X , p 2 o(f) = p 4 , and p 4 ° <j) = P2 



just using that the the function u a (9) introduced in (2.1 ) has period 7r. Furthermore, straightforward 
computations shows that 

cf> o n o (f> = T3 , (p o ts o (j) = n , cp o T2 o (j) = T4 , and c\> o T4 o </> = t 2 . 

Given an index k £ {1, .... 4}, we define its conjugate index k in the following manner: 

1 = 3, 3 = 1, 2 = 4, and 4 = 2. 



On the other hand, (6.1) is equivalent to 

g(x) = g(Tj {n(x))) if p, [x)pj (r, (a;)) > 0. 



(6.2) 



We have to check that this property holds for go (p. So let x £ [0, 7r] such that pi (x) pj (n (x)) > 0. 
Since Pi = p^o <j) and 4>° T -. — T i° 4> f° r an y * = 1) ■■■7 4, we have that Pi (x) pj (r^ (x)) > implies 



Pi (x) = pjo (p(x) > 

Pj (n (x)) = p~ } o 4> o n (x) = p~. o 7% o 4> (x) > 

Therefore, using the invariance of g expressed in (|6.2|) 



{Pj- (P] o (x) > 0. 



9(H X )) = 9 



Tj O Tr O 0(X 



so g o (j> also satisfies (6.2 ) and is S 2 -invariant. 



Proposition 7 Let g £ L 1 ([0, tt], fx) and suppose that 7rJ(<7) is invariant by S 2 and (f>. If a in (M2) 
and ( Ml ) is irrational, then g is constant. 

Proof. Let [0, tt]/ ' be the quotient space obtained by identifying ~ tt. This space is homeomor- 
phic to the unit circle S 1 . The function g : [0, tt] — > K induces a function [g] : [0, tt]/ ^^ — {[0]} — > M 
on the quotient that is well defined except at [0], the equivalence class of 0, because maybe 
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<?(0) 7^ 9(71")- We are going to prove that [g] is invariant by an irrational rotation. Since con- 
stants are the only integrable functions invariant by irrational rotations, [g] must be constant and 
g too. 

Recall that 713 (g) = tt^ (g) ° S 2 reads as 



g(x) = g{T j (T i (x))) i{p i (x)p j (n(x))>0. ^2) 

We are going split the interval [0, tt] into subintervals and, for each of them, we will explore the 
invariance properties of g implied by (6.2). 

(i) If x £ (0, tt — 4a), then pi (x)pi (ti (x)) = p\ (x)p\ (x + 2a) > so 

g (x) = g {t x (t x (x))) = g (x + 4a) . 

(ii) If x £ [tt — 4a, tt — 3a), then p\ (x)p2 (j\ (x)) — p\ (x) P2 (x + 2a) > 0. Writing x = tt — 4a + z 

with z £ [0, a), 

g (-7T — 4a + z) — g (j2 (ti(7t — 4a + z))) = g (tt — 2a — z) . 
Now, since g is (^-invariant, 

g (tt — 2a — z) = g (<f> [it — 2a — z)) = g (2a + z) 

Therefore, 

[g] ([x]) = [g] ([x + 6a]) if x <G [tt — 4a, tt — 3a). 

(iii) If x £ [tt — 3a, tt — 2a) then P2 (x) p% (T2 (x)) > 0. Writing x = tt — 2a — z with z £ [0, a), 

g (tt — 2a — z) = g (T3 (t2(tt — 2a — z))) = g (tt — 4a + z) . 

Using g = go(f), 

g (tt — 4a + z) = g ((f) (tt — 4a + z)) = g (4a — z) . 

That is, 

[g] ([ x ]) = [g] ([x + 6a]) if x £ [tt - 3a, tt - 2a). 

(iv) Suppose x £ [tt — 2a, tt — a), x = tt — 2a + z with z £ [0,a). Since P2 (x)p\ (T2 (x)) > on 
[7r-2a,7T-a),by Q, 

=g(r 2 oTi(a;)), r 2 o n(x) = 7r - z. 

Since g is invariant by (f>, 

g(n - z) = g (<p (tt - z)) = g (z) , z £ [0, a). 
Summarising, g (tt — 2a + z) — g(z) which implies 

[g] ([x]) = [g] ([x + 2a]) if x £ [tt - 2a, tt - a). 
Now z £ [0, a) and we can use the S 2 -invariance expressed in the first item so that 

[g}([x]) = [g}([x + 2a}) = [g]([x + 6a]). 
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(v) Finally, let x G [n — a, tt), x = n — a + z with z G [0, a). Since ps (x) P2 (ra (x)) > on [n — a, tt) 

g{x)=g (r 2 o r 3 (x)) , t 2 o t 3 (x) =ir — a-z. 



by (6.2). Using the 0-invariance of g, 

g(ir - a - z) = g ((/) (tt - a - z)) = g (z + a) . 

That is, g (x) = .g(z + a), which again implies 

[ff] (N) = M {[x + 2a]) if x € [tt - a,7r). 

Now x + 2a € [a, 2a) and 2a < tt — 4a as by assumption a < 7r/6; we are in the situation of 
(i) so 

[g}([x]) = [g}([x + 2a}) = [g}([x + 6a\). 

In conclusion, if g is invariant by S 2 and </), then 

ufui _ / Iff] ([i + 4a]) if x G (0, tt - 4a) , , 

[5JU Jj ~ I [ff]([x + 6a]) ifxG [7r-4a,7r). 

Let I; £ N such that (4k + 6) a > 7r > (4 (fe — 1) + 6) a. Observe that fc > 1 because a < 7r/6. 
Applying the invariance expressed in (6.3) iteratively, we have 

[g] ([x]) = [g] ([x + (4fc + 6) a]) for any x G (0, tt) . (6.4) 

In other words, [5] is invariant by the rotation of angle (4k + 6) a. If a/2ir is irrational, <? must be 
constant. ■ 



Proposition 8 Let C G B ([0, tt]) and suppose that 7r|(lc) is S 2 -invariant. If a in (M2) and (Ml) 
is irrational, then \x (C) equals either 1 or 0. 

Proof. Suppose that < (J,(C) < 1. We want to prove that ^(lc) cannot be ^-invariant. Since 
constants are trivially ^-invariant, we can subtract /i (C) from lc so that lc — /U (C) is still S 2 - 
invariant and has expectation 0. Let g := lc — A* (C) an d define g := | (g + g o 0), which is clearly 
0-invariant because 2 = Id. Since ^(g) is S'-invariant, so are 712(170 <j>) (Proposition [6| and 7r|(g). 
By Proposition [7j g is constant and equal to its expectation E M [g]. But [go^] = E M [g] = 
because dfi = sin (x) dx and sin (x) are invariant by <j>. Therefore, E^ \g] = and g, which is 
constant, must be equal to 0. That is, 

9 = -g° <t>- (6-5) 

Since 

1 - fi (C) > if x G C 



g(x) = l c (x)-^(C) 



-\i (C) < if x i C, 



(remember that we assumed < /x (C) < 1), we conclude from (6.5 ) that 1 — /x (C) = — (— [i (C)) 
M(C) = 1/2 and 

g{x) = 1/2<s=^xgC, 
g(x) = -1/2 «i6(7 c . 

Moreover, looking carefully at the proof of Proposition [7J we have 



[5] ([x + 4a]) if x G (0, 7r - 4a) 
— [g] ( [x + 6a] ) if x G [tt — 4a, 7r) , 
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where instead of g = g o <f> we have now used g = ~g o <f>. Therefore, (see (6.4)) 
[g\ ([«]) = - [9] (N + (4* + 6) a]) for any z G (0, vr), 



(6.6) 



where fc > 1 is such that (4fc + 6) a > 7r > (4 (A; — 1) + 6) ot. The minus sign in (6.6) tells us that 
the rotation Rp of angle (3 := (4fc + 6) a — ir sends C to C c and vice versa. 

Let lc = X)^L~oo a " e 2mx and lc= = e 2 " 11 the Fourier expansions of lc and lc c 

respectively. Since lc + lc c = l[o,7r]> we have that 

a n — —a n if n 7^ 



On the other hand, we already argued that lc Rp = lc c 
of lc- and lc?e are unique, we deduce that 



Imposing that the Fourier coefficients 



On e ^ flrj. 



-a n for any n ^ 0. 



As /3/2tt is irrational, e 2m ' 3 ^ -1 for any n^O, which implies that a n = for any n ^ so lc is 
constant a.s.. But this is clearly a contradiction because (6.5) implied /i (C) = 1/2. ■ 



7 Exactness of S 

In this section, we are going to prove that S is exact. This will imply Knudsen's strong law for the 
random billiard introduced in Section [2j 

Let T : fl — > fl be a measurable transformation of a probability space (ft, J 7 , p). From the 
measurability of T, we have the chain of c-algebras 

... C T~ n (J 7 ) C ... C T~ 2 (J 7 ) C T- 1 (J 7 ) C T. (7.1) 

The map T is called exact if the er-algebra Q := P|n>o (•^ r ) 0Iu y contains sets of measure either 
or 1. It is not difficult to prove that sets in Q are characterised by the property 

AeG ^ A = (T n {A)) for any n e N. 

This characterisation, in turn, implies that a map T is exact if and only if 

lim li{T n (A)) = 1 



for any i £ J such that p (A) > ([T31 Section 2.2]). Equation (7.1 ) reads at the level of L 2 spaces 
as 

... C U£L 2 C ... C U^L 2 C U T L 2 C L 2 

where L 2 = L 2 (f2, J 7 , jit) and t/r (/) = / ° T is the Koopman operator defined on L 2 . If a map 
is exact, then 

fl°° U£L 2 = C, 

which means that Uj, (/) — >• E M [/] in L 2 as n 00 (see [T31 Section 2.5]). If v is a probability 
absolutely continuous with respect to /1 with Radon-Nikodym derivative / € L 1 , then 

T> (A) = E M [/[/? (1 A )] -* fi (A) E M [/] - p(A), 
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that is, the sequence of measures T™v converges (weakly) to /i for any absolutely continuous measure 
v < n (d Proposition 4.4.1 (b)], [HI Section 2.6]). 

We want to show that the skew-type representation S : [0,1) x [0, 7r] — > [0,1) x [0,7r] of the 
model introduced in Section [2] is exact with respect to the Liouville measure to conclude that, for 
any absolutely continuous measure v <C (i on B ([0, tt]) and any A e S([0, 7r]), 

v {n) {A)=( l A (n 2 oS n )du= [ U~s (l -i, A ^\ d\®v — > n(A), (7.2) 

i.e., that the strong Knudsen's law holds. Observe from (7.2) that we only need to consider the 
action of the Koopman operator Us on square integrable functions that are the pull-back by TT2 
of functions in Lg ([0, tt], B ([0, tt]) , fx). We are going to denote L 2 , ([0, tt], B ([0, tt]) , //) simply by 
I? ([0, tt]) for the sake of a clearer notation. Consequently, we do not need to check that S is exact 
for the Borel cr-algebra of [0,1) x [0,tt] but it is enough to consider a smaller one, namely, the 
smallest cr-algebra that makes both the functions in ix* 2 (L 2 ([0,7r])) and S measurable. 

Definition 9 For any n £ N, let 

T n ^({^(lOfVn-n I CeB([0,7r])}) , 

where Ji x ...i n are as in Definition^ For any U € B ([0, tt]), the sets tt^ 1 (U) f] J^...^ and tt^ 1 (U) 
will be called generators (of ' J-„ ). 

The sequence of cr-algebras {.Fn} ngN define a filtration, 

JlCJ 2 C ... C T n C ... 

Indeed, for any U € B ([0, tt]) and any n > 1, a generator tt^ 1 (U) f] Ji 1 ...i„ can be written as 

N 

{u)f\j ix ,.. in = u (^(Eofy**™*-) ' 
fc=i 

which implies J-„ C J-n+i. Wc define T as the limit of this filtration. 
Definition 10 



Proposition 11 5 : [0, 1) x [0, tt] — > [0, 1) x [0, tt] is J 7 -measurable. For any f € L? ([0, 7r]), ^ (/) 
is aZso T -measurable. 

Proof. The second statement of the proposition is obvious. To prove that S is ^-measurable 
it is enough to show that S 1-1 (A) g T for the generators of T . So let A E T such that A = 
7T 2 1 (U) Pi J*i —in wnere U G B([0,7r]) and n > 1. Then, by definition of the sets Jj 1 ...i n 

AT 



fc=i 



which clearly belongs to J 7 . ■ 

The following two lemmas aim at getting a better insight of the structure of the cr-algebra J- '. 
As an immediate consequence of them, it will be enough to show that S is exact by looking at its 
action on the generators of T . 
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Lemma 12 Any set A € T n can be expressed as a finite union of disjoint generators. 

Proof. First to all, observe that, for a fixed n € N, the sets J^...^ are disjoint and form a partition 
of [0, 1) X [0,tt]. 

The complement of a generator tt^ 1 (U) |"| Ji 1 ...i„ is 

(V([/)fyn...*„) c = I y K2 1 (u)P\j jl ... jn \ ijK'tfMV), 

yi#ii,...,j«^i„ / 

where 1+J denotes the disjoint union, so it can be expressed as a finite union of generators. On the 
other hand, if we take a countable union of generators {C' / t} fc>1 , Ck = t^ 1 (Uk) H ^ , 

where Z^- ^ := {k > 1 | tt^ 1 (U k ) f| is in {Cm} m >i}- But 

where the indices i\, ...,i^ can only take finite number of possibilities. Therefore, a countable union 
of generators reduces to a finite union. ■ 

The union lj n >i JF n of a filtration J 7 ! C J 2 C ... is an algebra of sets. The c-algebra generated 
by an algebra of sets A is characterised by containing all the sets in A and the limit of all monotone 
sequence of sets. That is, a (A) coincides with the monotone class generated by A ([5J Section 1.3, 
Theorem 1]). This observation is the key to proving the following Lemma: 



Lemma 13 If Ae T =a ^Un>i-^nJ ^ as P os tti ve probability, then A contains a generator of posi- 
tive probability. 

Proof. Let A g cr (Uri>i^ r ™)- If ^4 € Un>i^"«> tnen ^ e -F n f° r s °me n e N. By Proposition 



12 A can be expressed as a finite union of generators and at least one of them must have positive 
probability. If A ^ |J n>1 J n , then there exists a monotone sequence of sets {B n } nl - K C lj n >i Z n 
such that A = lim„_B„. We are going to deal with the cases {B n } neN increasing or decreasing 
separately. 

• If {B n } neN is increasing, then A = U„>i B n . Since A ® ji (A) > 0, at least one B n must have 
positive measure, A ® /x (-B„) > 0. But _B„ g ,F m for some m, so it contains a generator of 
positive probability. 

• Let {B n } neN C U n >i-^n be a decreasing sequence such that A = f] n> i B n . Suppose that 
< \l (A) < 1 (otherwise there is nothing to prove). We have 



A =(u„> 1 s ») c =((y n>1 ^ + A^)y^) c , 
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so that A can be expressed as the complement of a countable union of disjoint sets. Rename 
d := B\ and C n := B c n+1 \B c n , and express any C n as a disjoint union of generators (Lemma 



121 



In this decomposition, we only consider sets of strictly positive probability. Then, 

I TTn 1 (U r _) H J« i» = l+l (l I 



n>l r n = l 



^ c = i+i i+i ^2 1 (^j n ^ ».*? = i+j u me/ ,„ .„ ^ (t^m) n j <? -*? ( ? - 3 ) 



where := {m > 1 | 7r 2 1 (J7 TO ) p| J»»...i» C yl c = |+| n C„}. We claim that there exists a 

generator G of strictly positive probability that does not intersect A c , i.e., G C. A. 



Let be a sequence appearing in the decomposition (7.3 1. If A*(Um£7n „ < 1) then 

is a generator of strictly positive measure that does not intersect A c . If MUme/n , n U m ) = 1 



for any finite sequence in (7.3 1, then A c — (+) i „ in Ji«...^ a.s. and since we assumed 

that < /i (A c ) < 1, there must exist some set Ji 1 ...i k with positive probability that does not 
appear in the decomposition of A c . That is, Ji t ...i k — tt^ 1 ([0, 7r]) p| Ji x ...i k C A 



Theorem 14 for any A G J 7 , 

lim A®ju(S n (A)) = 1. 

n— >oo 

In other words, S : [0, 1) x [0, 7r] — > [0, 1) x [0, 7r] is an exact endomorphism o/([0, 1) x [0, 7r], J 7 , A <g> /i). 



Proof. Let 4 e J be such that A ® fj, (A) > 0. By Lemma 13 A contains a generator G = 
7T 2 1 (U) Pi Jii—i„ of positive measure for some U £ B ([0, 7r]), n > 0. It is an immediate consequence 
of the definitions that, after n iterations, the set 5™ (G) has all its fibers of length 1, that is, 

A (([0, 1) x {x}) P| S n {G)) = 1 for any x = r il o ... o r in (it) , it £ [/. 



Let B := S" (G) and let I x := [0, 1) x {x} be the fiber at x e [0, tt]. Looking carefully at (Ml ) 



and (M2), one can see that 



t 2 o t 2 — t 4 o t 4 — T\ o r 3 = r 3 o n = Id, 

where Id denotes the identity on [0,7r], and 

p 2 (x)>0^ p 2 (x) p 2 (t 2 (x)) > 
Pi (x) > => p 4 (#) P4 (i"4 (a;)) > 
pi (a:) > =>> pi (x)p 3 (n (at)) > 
p 3 (x)>0^ p 3 (x)pi (t 3 (x)) > 0. 
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These remarks imply that, if I x l~l Ji ^ 0, 

i x n J t c s 2 (/ x n Ji) . 

Since all the fibers of B have length 1, B C S 2 (B) and, iteratively, 

BCS 2 (B) C ... C S* 2 " (B) C ... 

Define C := U„> S 2 " (#)• Obviously S 2 (C) = C, the fibers of C have all length 1, and A®/i (C) = 
lim n _ ) . 00 A ® (J, {S 2n (B)) . We want to show that A <g> /i (C) = 1. Observe that A®/j (C) > because 
A <g> /i (B) > 0. From S 2 (C) = C we have 

C C S" 2 (S 2 (C)) = S- 2 (C) ; 

but C and S 1-2 (S 2 (C)J have the same measure (S 2 is A ® /i-preserving), so lc is a S^-invariant 
function a.s.. Since lc = AsK^-WC)) a - s - (the fibers of C have all length 1 a.s.), 7T2 (C) has full 
measure by Proposition [8j Therefore, A ® /z (C) = 1 and 

lim S 2 " (-B) = 1. 

n— >oo 

This proves that S* 2 is exact. In general, we have 

S 2n (B) C 5- 1 (5 (S 2 ™ (B))) C S- 1 (S 2n+1 (B)) 
S 2n+1 (B) C 5- 1 (5 (S 2n+1 (B))) C 5- 1 (s 2 <" +1 > (B)) 

so that, using again that S 1 is A ® /i-preserving, 

A ® ^ (5 2n (B)) < A ® ^ (S 2n+1 (B)) < A ® ^ (s 2 ^ 1 ) (B)) , 

which proves that {A <£> fi (S n (B))} ngN is increasing and converges to 1, i.e., S is also exact. ■ 

To sum u p th e ske w- typ e representation S : [0,1) x [0, n ] — > [0,1) x [0, w] associated to the 
random map (Ml) and (M2) is exact which implies, by (7.2 1, that u^ n > — > fi as n — > oo for any 
initial distribution v on [0, it] absolutely continuous with respect to the law fi (A) = \ j A sin (9) d9, 
A G B ([0, 7r]). Since sin (9) > on (0, ir), fi and the Lebesgue measure A are absolutely continuous 
with respect to each other on (0, n). In other words, strong Knudsen's law v^ n ' — > \x holds for any 
initial distribution absolutely continuous with respect to A. 



8 Knudsen's law. Simulations 



In this section, we are going to show that numerical simulations are in accordance with theoretical 
results. With this aim, we take an arbitrary initial distribution on [0, 7r] such as in the following 
picture and we make it evolve according to our random system (Ml) and (M2). 
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D.5 1 1.5 2 2.5 3 

Angle (45 bins) 

Initial normalised distribution of particles. 



Since we can only simulate a finite number particles in a computer, in this experiment we take a total 
amount of 30000 balls and divide [0, 7r] in subintcrvals of the same length (45 in our experiment). 
That is, we approximate the initial distribution by a step function. In each subinterval, we put 
the proportion of balls according to the probability density function above. The initial angle 
associated to any of those balls is the middle point of the interval where they fall. The simulation 
then goes as follows. At the i-th step, we take the j-th particle with angle Qy' and a random 
number uniformly distributed between and 1, one for each particle. If (yf\o^) £ then 

Of = TkiOj^), where {r fc } fc= i,...,4 are as in 
the distribution stabilises as follows: 



0.5 
045 

0.4 
035 

0.3 
025 

2 
15 

0.1 
005 


0.5 1 1.5 2 2.5 3 

Angle (45 bins) 

Final normalised distribution of particles. In red, \ sin(x). 

As we can see, the outline of the final distribution tends to the graph of |sin(:r). The small 
inaccuracy is explained by the fact that only a finite (i.e., a discrete) number of initial angles is 



Ml), and so on. After a long number of iterations, 
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considered. The smaller the subintervals in which [0, n] is divided are, the better the final distri- 
bution approximates \ sin (a;). We have repeated this experiment over several initial distributions 
obtaining always similar results, which experimentally confirms the validity of Knudsen's law for 
our model. 



A Appendix 

Proof of Lemma [3} To start with, will prove that 

S Ui P| J^-.-i^j = I Tin {x) P| Jii---i n -i ■ 
From the very definition of S and J^...^, we have 

S Ui P Jii - inJ ^ In n (x) P Jit—in-! ■ 

On the other hand, let uj £ I Ti (x) D ^ii— £«-i- Then u> — (y,Ti n (x)) for some y £ [0,1). Let 
lj' = (pi n (x)y + Y^jLi Pj( x )i x). It is obvious that u>' <E Ji n as pi n (x) > because I x f] Ji 1 ..-i n is not 
empty Moreover, S (uj') = lj. Therefore to £ S (I x f] Ji x —i n ) and 

5* Ui P Jii— i„J 2 ^n n {x) P Jh—i n --L ■ 

Identifying both fibers I x |"| Ji 1 ---i n and I Tin (x) fl Jh—i„-i as subsets of [0, 1), the map 

y ' — > v>in (y> n n (x)) = ^-j^ (y - Y?£iPj ( x i) > 

is then a diffeomorphism. After having identified I x f] Ji 1 ...i n as a subset of [0, 1), we can compute 
its Lebesgue measure A (I x f] Jii--.j n )- Consequently, 



.(/*fy<i" 



/ p ln (x)d\ = p in (x) j dX 



= ft„(!c)A(L )|iW f|Jj 1 ..Vi) • (A.2) 
Applying iteratively ( |A.2[ ), we obtain 

A (lxf]Jii-ir^ =Pt„(»)Pt„_i ( r i„(a;)) ■ • - Pi faa o • • • o Ti„(a;)). 
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